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ABSTRACT 


The  location  of  three-dimensional  separation  on  bodies  such  as  submarines  is  im¬ 
portant  diagnostic  information.  Considerable  changes  in  the  side  forces  and  moments 
occur  with  the  undesirable  phenomenon  of  separation.  Three-dimensional  flow  separations 
occur  very  close  to  local  minima  in  the  magnitude  of  the  wall  shear  stress.  A  1000-flush- 
surface-sensor  constant-current-anemometer  separation  location  system  is  described  here 
for  locating  these  local  minima  along  the  axis  of  a  body. 

As  shown  here,  the  convective  heat  transfer  coefficient  from  a  sensor  is  proportional 
to  the  one-third  power  of  the  local  skin  friction  magnitude  for  the  directionally-insensitive 
sensors  described  here.  The  constant-current  anemometer  design  permits  the  use  of  strips 
of  serially-connected  sensors  that  require  only  one  current  source  per  strip.  These  compact 
electronics  can  be  housed  inside  the  body. 

The  voltage  drop  across  each  sensor  is  measured  for  a  number  of  different  computer- 
controlled  c\u:rents  id  related  to  the  local  heat  transfer  coefficient  and  the  fixed  parame¬ 
ters  for  that  sensor.  Each  sensor  mounted  on  the  body  forms  a  first  order  dynamic  thermal 
system  whose  time  constants  can  be  experimentally  determined.  Thereafter,  such  as  in 
transient  flow  separation  conditions,  measurements  only  for  one  current  are  necessary  for 
a  given  pre-calibrated  sensor. 

Results  are  presented  for  two  sets  of  tests:  (1)  flow  in  the  VPI&SU  Water  Tunnel  over 
a  20°  angle-of-attack  circular  cylinder  and  (2)  a  6:1  prolate  spheroid  at  15°  angle  of  attack 
towed  in  the  David  Taylor  Research  Center  140’  Basin.  The  separation  location  results 
are  in  close  agreement  with  those  from  earlier  wind  tunnel  tests  by  other  investigators. 


I.  Introduction 

In  general,  premature  flow  separation  is  an  unwanted  phenomenon  because  it  leads 
to  less  than  ideal  performance  of  a  flow  machine.  In  the  case  of  a  submarine,  crossflow 
separations  form  on  the  leeside  of  the  stern  during  a  turning  manuever.  Large  streamwise 
vortex  structures  that  contain  much  circulation  are  a  result,  causing  a  considerable  change 
of  the  side  force  and  the  vertical  turning  moment  of  the  submarine.  The  behavior  during 
dynamically  unsteady  manuevers  is  even  more  compUcated  because  of  the  time-dependent 
locations  of  the  separations  and  strengths  of  the  shed  vortex  structures. 

By  ’’separation”,  we  mean  the  entire  process  of  ’’departure”  or  ’’breakaway”,  or  the 
breakdown  of  the  boundary-layer  flow  concept.  An  abrupt  thickening  of  the  rotational-flow 
region  next  to  a  wall  and  significant  values  of  the  normal-to-wall  velocity  component  must 
accompany  breakaway,  or  otherwise  this  region  would  not  have  any  significant  interaction 
with  the  free-stream  flow.  In  three-dimensional  separations  the  rotational  viscous  layer 
can  depart  from  the  wall  without  the  surface  shear  stress  necessarily  falling  to  zero,  such 
as  occurs  in  two-dimensional  separations.  In  three-dimensional  separations  the  wall  shear 
stress  is  zero  only  at  the  singular  points  (Simpson,  1985). 

Figure  1  shows  a  three-dimensional  open  sep>aration  with  several  characteristic  fea¬ 
tures.  Limiting  streamlines  rising  on  either  side  of  the  line  of  separation  are  prevented 
from  crossing  the  presence  of  a  stream  surface  stemming  from  the  line  of  separation 
itself.  The  surface  skin  friction  lines  on  each  side  of  a  separation  Hne  converge  toward 
one  another,  indicating  that  the  low-velocity  near  wall  viscous  flow  on  each  side  of  a  sep¬ 
aration  line  accumulates  along  that  line  before  moving  away  from  the  wall.  There  is  no 
flow  reversal  along  the  separation  line.  The  skin  friction  does  not  approach  zero,  although 
it  reaches  a  local  minimum  where  the  lowest  velocity  near  wall  fluid  is  located  (Tobak 
and  Peake,  1982).  Appendix  A  is  a  more  detailed  discussion  of  the  relationships  between 
surface  parameters  and  the  separation  line  location. 

Figure  2  shows  data  from  wind  tunnel  tests  of  a  6:1  prolate  spheroid  that  indicate 
two  separation  locations,  SI  and  S2,  approximately  where  local  minima  of  the  skin  friction 
magnitude  occur.  Note  that  the  skin  friction  lines  converge  to  the  separation  line  from  ei¬ 
ther  side,  as  previously  discussed  above.  Therefore,  a  technique  which  is  to  experimentally 
locate  three-dimensional  separation  should  either  find  the  location  of  the  minimum  surface 
shear  stress  magnitude  as  discussed  in  Appendix  A  or  give  the  direction  of  the  shear  stress 
lines  so  that  their  convergence  can  be  identified. 

Moskovitz  and  Simpson  (1990)  investigated  several  schemes  for  detecting  one  of  these 
two  features.  Since  the  locations  of  the  separations  change  with  angle  of  crossflow,  speed 
of  the  submarine,  and  the  unsteady  effects  of  a  particular  manuever,  a  separation  loca¬ 
tion  system  should  have  suflScient  transient  response  to  capture  the  changing  locations  of 
separation.  They  concluded  that  a  series  of  flush-mounted  heated-foil  elements  was  the 
best  approach  for  locating  minima  in  the  skin-friction  magnitude  in  time-dependent  water 
flows. 


This  report  describes  some  of  the  effort  to  date  in  the  development  of  a  constant 
current  anemometer  system  to  determine  the  locations  of  three-dimensional  crossflow  sep¬ 
arations  on  submarines.  The  following  section  is  an  analysis  that  shows  the  relationship 
between  forced  convection  heat  transfer  and  the  local  mean  skin  friction.  Section  III  of 
this  report  is  devoted  to  natural  or  free  convection  heat  transfer  effects  that  can  become 
important  at  sufficiently  low  flow  velocities. 

A  specially-designed  heat  transfer  element  that  is  insensitive  to  flow  direction  is  de¬ 
scribed  in  section  IV,  along  with  data  reduction  equations  and  an  analysis  of  the  transient 
response  of  a  sensor. 

Section  V  deals  with  sensor  covering  and  waterproofing  materials. 

Section  VI  describes  the  electronics  that  have  been  designed  to  accommodate  up  to 
1000  sensors  that  will  be  mounted  on  the  stern  of  the  16’  long  rotating  arm  model  for  use 
at  the  David  Taylor  Research  Center.  Portions  of  these  electronics  were  modified  in  order 
to  conduct  the  tests  in  the  VPI&SU  Water  Tunnel  and  in  the  DTRC  140’  Carriage  Basin 
that  are  described  here. 

Section  VII  describes  the  experimental  results  of  tests  performed  in  the  VPI&SU 
Water  Tunnel.  A  series  of  tests  was  conducted  on  three-dimensional  flow  over  a  circular 
cylinder  at  22  degrees  angle  of  attack  to  the  flow.  Section  VIII  describes  a  series  of  tests 
conducted  at  the  David  Taylor  Research  Center  on  three-dimensional  flow  over  a  6:1  prolate 
spheroid  at  15°  angle  of  attack.  These  tests  show  that  data  from  this  constant-current 
anemometer  system  locate  the  of  3-D  separations. 


II.  Analysis  of  Forced  Convection  from  a  Single  Heated  Sensor  to  a  High  Prandtl  Number 
Fluid 

For  a  small  surface  area  constant  current  anemometer  sensor  used  in  water,  we  can 
employ  a  high  Prandtl  number  analysis  that  restricts  the  thermal  boundary  layer  to  the 
viscous  sublayer  next  to  the  wall.  The  approach  velocity  profile  is  already  developed  and 
does  not  change  much  over  the  length  of  the  sensor.  Here  we  wish  to  obtain  the  relationship 
between  the  mean  velocity  gradient  at  the  wall  and  the  average  heat  transfer  rate  over  the 
sensor. 

The  thermal  energy  transport  equation  near  the  wall  is  given  by 


U 


dx 


(1) 


for  locally  two-dimensional  flow.  Since  the  flow  upstream  of  the  sensor  is  unheated  the  in¬ 
coming  flow  temperature  is  the  same  as  that  of  the  freestream.  Too  •  The  thermal  boundary 
layer  grows  from  the  upstream  edge  of  the  sensor,  but  remains  entirely  within  the  near  wall 
region  if  the  sensor  is  short  enough  in  the  streamwise  direction.  We  make  that  assumption 


here  which  justifies  the  use  of  the  following  near  wall  velocity  profile 


U 


dy) 


y 


(2) 


where  dU/dy)w  is  the  mean  velocity  gradient  at  the  wall.  The  surface  shear  stress  is 
given  by  Ty,  =  fj,{dU/dy)y,. 

Substituting  the  non-dimensional  temperature 


T  = 


T-Ty, 
Too  ~  Tu, 


(3) 


into  equation  (1)  and  using  (2)  produces 


dU\  dT  d^T 

9y)y,^dx  ~  ^dy^ 

(4) 

The  boundary  conditions  on  the  fluid  for  this  problem  are 

T  =  0  at  y  =  0 

(5a) 

T  =  1  as  y  — »  00 

(56) 

T  =  1  at  X  =  0 

(5c) 

A  dimensionless  similarity  variable  rj  can  be  found  that  reduces  equation  (4)  to  a 
linear  ordinary  differential  equation 

v(du/av)i'^ 

^  (3ax)^/^ 

(6) 

T"  -1-  rf^V  =  0 

(7) 

Using  equation  (6)  the  boundary  conditions  then  become 

H 

II 

o 

-3 

II 

O 

(8a) 

T  =  1  as  77  oo 

(86) 

The  solution  to  equation  (7)  using  boundary  conditions  (8)  is 

T  =  M  r 

Jo  1 

(9) 

where 

1/M  =  j  =  |r  =  0.893 

(10) 
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The  heat  flux  at  the  surface  is  given  by 


)  =  -KT^ 

/  tu 


(11) 


which  becomes 


using  equation  (9). 


=  Mk{T^  -  Too) 


(du/dyfj^ 


(12) 


For  a  rectangular  sensor  of  constant  width  w  in  the  spanwise  direction,  we  can  obtain 
an  average  heat  transfer  rate  over  its  length  £ 


^Mkw(T^  -  r„) 


(dU/dy)]!^ 

(3a)>/3 


which  means  that  the  average  convective  conductance  is 


h  = 


Q 

ew{T^  -  Too) 


{au/ay)U^ 


(13) 


(14) 


Note  that  the  measured  h  is  directly  proportional  to  (dUfdy)^^^  and  all  other  parameters 
in  this  equation  are  properties  of  the  fluid  and  the  sensor. 

Thus,  if  the  sensor  is  sensitive  only  to  the  magnitude  of  the  wall  shear  stress,  as  in 
the  present  design,  measurements  of  h  from  a  row  of  nearby  sensors  permit  the  location 
of  local  minima  in  h  and  {dU/dy)w  Here  we  calculate  ^  and  {T^  -  Too)  directly  from 
electrical  measurements  of  the  sensor  condition,  as  described  in  the  section  below  on  the 
sensor  design  and  behavior. 

III.  Consideration  of  the  Nature  of  Combined  Forced  and  Free  Convection  Heat  Transfer 


The  analysis  in  the  last  section  assumed  that  free-convection-induced  heat  transfer 
was  negligible.  Clearly  this  is  not  always  the  case,  especially  near  stagnation  points  and 
zones  of  two-dimensional  separation  where  there  is  flow  reversal.  In  this  section  an  integral 
analysis  of  combined  forced  and  free  convection  will  be  presented  that  will  allow  us  to  assess 
the  relative  effects  of  free  convection  for  three-dimensional  separations. 

The  governing  thermal  energy  equation  in  the  fluid  is 


dx  dy  dz  dy^ 


(15) 


We  will  allow  forced  convection  in  the  x  direction  and  the  positive  z  direction  opposite  to 
the  gravitational  acceleration  vector  g.  As  in  the  last  section,  we  will  assume  that  equation 


(2)  holds  and  V  v  0  near  the  wall.  We  will  utilize  the  commonly  used  free-convection- 
induced  velocity  profile  (Holman,  1990) 


where 


Here  P 


1  dp 
~p  W 


6t  is  the  local  thermal  boundary  layer  thickness,  and  i/  is  the 


kinematic  viscosity.  Equation  (16)  satisfies  the  exact  boundary  condition  for  d'^Wfdy'^  at 
the  wall  from  the  z  direction  momentum  equation,  the  no  sUp  condition  W  =  0  at  t/  =  0, 
and  the  artificial  condition  dWfdy  =  0  at  y  ~  St- 


We  will  use  the  temperature  profile 


T-T^ 


(18) 


Tw  '^oo 

which  satisfies  the  conditions:  T  =  at  y  =  0;  and  T  —  Too  and  dT/dy  =  0  at  y  =  St- 

Substituting  all  of  these  relations  into  equation  (15)  and  integrating  with  respect  to 
y  produces 


pcp  ^  {T^  Too)S^ 


,d 


/.'(■-i)  Hi), 


+p  Cp  G  St  (Tu,  —  Too)  ^ 


Of-  .  .1  « 

—  ~  Qw 

Ot 

for  a  constant  surface  temperature  sensor.  This  further  becomes 


1. 

12  dy  dx' 


■  ^  9^w  Too)  _^/c3'v  _  ^ 

10  a*.  ion  ..  ^ 


120  u 


— 

St 


After  further  rearranging,  this  equation  becomes 
7  pCp  dU  dS-^  p  0  Cp  g 


dS^ 

36  k  By  dx  24o7T  "af  ”  * 


(19) 


(20) 


(21) 


The  initial  conditions  for  this  non-linear  differential  equation  for  a  rectangular  sensor 
of  length  i  and  width  w  are  St  =  0  along  x  =  0  and  St  =  0  along  2  =  0.  The  latter 


condition  is  not  exactly  rigorous  because  W  is  not  zero  along  2  =  0  because  of  continuity 
requirements  on  the  free-convection-induced  flow. 

Integrating  equation  (20)  over  the  area  of  a  sensor  produces  an  equation  for  the  total 
heat  transfer  rate 


7pcp  .du 


■  rW 

L  ' 


6t  ^dz 


J  x^t 


+  - 1207 -  X 


q  (22) 


with  the  first  term  representing  the  forced  convection  heat  transfer  and  the  second  term 
the  free  convection. 

The  property  groupings  for  these  two  terms  produce  insights  on  the  relative  magni¬ 
tude  and  temperature  dependence  of  forced  and  free  convection.  Over  the  range  of  16° C 
to  38°C, where  we  expect  to  operate,  p  and  Cp  for  water  are  constant  within  ±1/4%  while 
k  increases  less  than  8%.  Thus  the  thermal  properties  grouping  associated  with  the  forced 
convection  term  of  equation  (21)  decreases  about  8%  from  10°C  to  38°C.  The  free  convec¬ 
tion  grouping  p  ^  Cp  g  f  k  v  varies  almost  linearly  with  temperature  and  increases  by  a 
factor  of  about  5.2  over  this  temperature  range.  Thus,  the  coefficient  of  the  second  term 
of  equation  (22)  varies  with  about  (T„,  —  Tqo)^  while  the  forced  convection  term  coefficient 
varies  with  about  (T^,  —  Too)- 

The  closed  form  solution  of  equation  (21)  is  implicit  in  terms  of  br  and  is  not  readily 
available  to  substitute  into  equation  (22).  Since  we  wish  to  obtain  the  forced  convection 
contribution  to  equation  (22),  we  can  develop  an  experimental  procedure  using  data  at 
various  (T^,  -  Too)  and  the  form  of  equation  (22) 


Equation  (22)  can  be  written  in  the  form 

Fx  (T„  -  r„)  +  F,  (r„  -  T^f  =  5 


(23) 


Fi  is  the  forced  convection  term  coeflacient  that  varies  principally  with  dU/dy  since  the 
properties  do  not  vary  much  with  temperature  and  the  geometry  of  a  given  sensor  is  fixed. 
F2  is  a  function  associated  with  the  free-conwction  and  depends  on  the  fluid  properties, 
sensor  geometry,  and  sensor  operation.  It  increases  significantly  with  the  average  temper¬ 
ature  Tm  =  (Tu,  ±  Too)/2,  which  is  used  to  evaluate  properties.  Fi  is  the  slope  of  a  plot  of 
q  vs.  (Tu,  -  Too)  as  (T„  -  Too)  -»  0,  as  shown  in  equation  (23). 

In  order  to  estimate  when  free  convection  effects  may  be  significant,  we  can  look  at 
the  ratio  of  the  free  convection  term  to  forced  convection  term  in  equation  (21).  The 
dimensionless  ratio  of  the  coefficients  to  the  derivatives  is 


„  ^  (T^-Too)gr^  ,(J_P£p^ 

\  240  u  k  )  \36  k  dy 


(24) 


This  indicates  that  free  convection  is  small  when  6t  is  small  in  the  presence  of  any  forced 
convection.  We  can  estimate  6t  from  the  pure  forced  convection  solution  given  in  the  last 
section. 


6 


(25) 


Forced  convection  Sr 


2  (3  a 

M  (dU/dyy/^ 


0.0556 

This  indicates  that  we  can  minimize  free  convection  effects 
short  sensors  {£  is  small). 


Thus 


R  = 


/(dU/dyff^  (26) 

by  using  low  (T^  —  Too)  and 


Consider  water  at  the  average  properties  temperature  of  about  21°C,  then  equation 
(26)  produces 


(  0.63  \ 

°C  sec^/^j 


(du/dyyr-i 


(27) 


For  /2  <  0.1,  then  this  indicates  that 


dy 


> 


3.98 

mi/4  oca/-* 


(T„  -  e'!' 


(28) 


As  long  as  this  criterion  is  met,  free  convection  effects  should  be  less  than  about  10%  of 
the  total  heat  transfer. 


rV  Sensor  Design  and  Behavior 

figures  3  and  4  show  two  designs  of  nominally  directionally  insensitive  constant  cur¬ 
rent  anemometer  sensor  foil.  The  sensor  shown  in  Figure  3  have  a  ±  6%  variation  of  heat 
transfer  with  direction,  while  the  second  design  has  a  ±  1%  variation.  These  variations 
were  deterriilned  for  various  orientations  of  a  sensor  by  assuming  that  the  convective  heat 
transfer  from  a  given  length  of  the  sensor  element  varied  with  the  cosine  of  the  angle  of 
the  flow  from  the  normal  to  the  element;  the  contributions  from  each  element  were  then 
summed. 

The  sensor  foil  of  thickness  t  is  mounted  on  a  substrate  material  of  thickness  tg  and 
covered  with  a  material  of  thickness  tc,  and  has  a  planform  surface  area  A*.  A  layer  of 
adhesive  of  thickness  ta  is  between  each  of  the  other  layers.  Thus,  a  sensor  is  a  five  layer 
sandwich  as  shown  in  Figure  5.  When  applied  to  the  surface  of  a  model,  an  additional 
layer  of  adhesive  must  be  used.  For  the  MINCO  sensors  shown  in  Figure  3,  tc  =  tg  = 
ta  =  0.0005  inches  and  t  =  0.0002  inches.  The  cover  and  substrate  materials  are  kapton 
and  the  foil  is  Balco  (70%  nickel,  30%  iron)  with  the  temperature  coefficient  of  resistivity 
a  =  0.0051/0°.  Table  1  is  a  list  of  properties  of  materials. 

Sensors  are  manufactured  in  rows  with  interconnecting  leads  as  shown  in  Figure  6.  A 
constant  current  I  from  a  single  source  is  passed  through  the  entire  row  of  sensors.  The 
Joule  heating  in  each  sensor  I^R  is  convected  away  by  the  fluid.  Ideally,  we  would  prefer 
no  conduction  into  the  substrate  and  model  and  no  thermal  resistance  to  heat  transfer  into 
the  water. 


In  steady  state,  the  Joule  heating  equals  the  heat  transfer  into  the  substrate  and 
model  and  the  heat  transfer  into  the  water 


q^I^R  =  I^R^[l  +  a{Tf-T^)]  =  {Tf  -  T^)im  +  7^2)  (29) 


where  T f  and  Too  are  the  temperatures  of  the  foil  and  water,  respectively. 
Here, 

TTi  = 


As 


(fc  +  fe  +  ^) 

is  the  thermal  conductance  into  the  water  and 


7r2 


As 


ix.  -L  2ia.  -L  Im. 

1.  'V 

s  ftfri 


(30) 


(31) 


is  the  conductance  into  the  model  of  wall  thickness  tm  and  thermal  conductivity  k^.  For 
a  fiberglass  model  wall  of  the  order  of  1/4  inches  thick,  7r2  «  tti  x  10“^,  and  the  heat 
transfer  is  mainly  into  the  water.  Ideally,  the  thermal  resistance  of  the  cover  and  adhesive 
material  {tc/kc  +  ta/ka)  should  be  made  close  to  zero. 

Equation  (29)  can  be  rearranged  to  give 


(Tf  -  Too) 


(tti  +  7r2  -  aPRoo) 


(32) 


Because  R  varies  with  a  and  (T/  —  Too)  as  implied  in  equation  (29),  the  steady  state  voltage 
drop  across  a  sensor  is 

V  =  /flocIl  +  a(T/-Too)]  (33) 


Equations  (32)  and  (33)  can  be  combined  and  the  voltage  sensitivity  to  h  calculated 

dh  \dni)\ah)  -  aPR^yh'‘  '  ’ 

In  the  limit  of  7r2  — ♦  0  (no  heat  transfer  to  the  substrate)  and  xi  -+  hAs  (no  cover  and 
adhesive  thermaJ  resistance),  equation  (34)  becomes 


W  ^  -aPRl,  As  ,  , 

dh  (hAs-aP  R^)^  '  ’ 

In  both  equations  (34)  and  (35),  we  see  that  the  sensitivity  of  the  voltage  with  respect  to  the 
convective  conductance /i  varies  stronger  than  the  individual  variations  of  Q,  P,  R^,  ^/h^- 
and  I /As-  Equation  (35)  is  an  important  design  equation,  indicating  that  we  should  min¬ 
imize  heat  transfer  to  the  substrate  and  model,  reduce  thermal  resistance  of  the  covering 
materials  within  waterproofing  constraints,  select  a  high  a  foil  material,  minimize  As 


within  manufacturing  constraints,  and  maximize  I  and  Roo  within  maximum  voltage  drop 
and  manufacturing  constraints. 


In  practice  we  measure  the  cold  resistance  Roo  of  each  sensor,  the  voltage  V  across 
each  sensor,  and  the  current  /  which  is  held  constant.  Thus  we  can  calculate  the  foil 
temperature  difference 


and  the  heat  transfer  rate 


q  =  VI 


(37) 


Using  the  defining  equation  for  h 


q  =  h  -  Too) 


(38) 


in  terms  of  the  sensor  cover  outer  temperature,  assuming  that  all  of  the  heat  transfer  is 
into  the  water,  and  using  the  thermal  resistance  ratio 

T^-T^  =  _ I _ 

+  fe)  +  1 

we  can  calculate  the  surface  temperature  difference 


T„-T^  =  (T,  -  I"™)  - 


(40) 


Note  that  q  and  (T/  —  Too)  are  derived  from  experimentally  measured  quantities  and 
equations  (36)  and  (37)  and  the  other  quantities  axe  properties  of  the  sensor  and  the 
covering  materials.  Finally,  equation  (38)  can  be  used  to  determine  h. 


To  examine  the  transient  response  of  a  sensor,  we  can  consider  a  time-dependent 
contribution  to  the  convective  conductance,  ^(t),  which  forces  a  transient  temperature 
difference  d(r)  of  the  foil.  The  governing  linearized  first  term  unsteady  heat  transfer 
equation  is 


M  ^  =  -[tti  -H  7r2  -  qI^Roo]0  - 
dr 


(Tf-Too)n^hi(T) 


(41a) 


or 


M  dO/dr  =  —N  $  —  P  /ii(t) 


(416) 


where  tti,  7r2,  (T/  —  Too),  and  h  are  from  the  steady  state  mean  flow  case  described  above, 
and 

M  =  As  Y^(pi  CjH  ti  Bi)  (42) 

t 

is  the  sum  of  the  lumped  parameter  thermal  capacitance  for  each  layer  i  of  the  sensor 
sandwich.  Here  Bj  is  the  steady  state  ratio  of  (Tj  —  Too)/(Tf  —  Too)  for  each  layer  of  the 
sandwich. 


For  a  sinusoidal  variation  of  hi  (r) 


hiir)  =  h'  sin(a;T) 


the  solution  of  equation  (41)  is 


0(r)  = 


-  uPh' 


!JL  • 


sin(a;T)  —  cos(u;t)  +  e 


for  0  =  0  at  T  =  0.  For  large  time  r 


eiT)  = 


P  h'  sin(u;r  +  0) 
M(a;2  +  isr2/M2)i/2 


where 


0  =  TT  +  tan' 


.1  /  u)M\ 

[~ir) 


is  phase  shift.  Note  that  as  M  — ►  0,  9  is  180°  out  of  phase  with  hi.  The  voltage  variation 
across  the  sensor  is 

V(t)  =  alR^eir)  (47) 

For  a  ramp  variation  of  hi  (r) 

Hi{r)  =  (dh/dr)ir  (48) 

the  solution  of  equation  (41)  with  0  =  OatT  =  Ois 

For  large  time  r,  the  time  delay  D  is  found  by  equating  0{t)  for  Af  =  0  to  6{t  +  D)  for 
M^O: 

=  -p  (^)  _  X  =  _£  (*)  ^  (.  + 1,  _  I)  (50) 


Thus,  the  time  delay 


D  =  M/N 


which  indicates  that  piecewise  continuous  ramp  variations  of  hi  (r)  result  in  a  time  delay 
response  that  is  independent  of  the  ramp  slope  (dhi/dr). 


V  Sensor  Covering  and  Waterproofing  Materials 

The  analysis  in  section  IV  describes  the  need  for  a  thin  coating  with  high  thermal 
conductivity  to  give  the  sensor  the  highest  sensitivity  obtainable  for  a  given  design.  The 
coating  used  and  its  application  are  of  primary  concern  to  the  success  of  the  system.  For 
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this  reason  it  has  been  decided  to  manufacture  future  sensor  strips  without  a  cover,  since 
the  thinnest  kapton  or  upilex  covers  available  are  0.5  mils  with  a  thermal  conductivity  k  no 
better  than  1.4  BTU-in/(R-hrft^),  such  as  used  in  the  tests  reported  here.  Removal  of  the 
cover  also  eliminates  a  layer  of  adhesive  (on  the  order  of  0.5-1. 0  mils  with  a  conductivity 
of  1.3  BTU-in/(R-hr-ft^)).  In  the  future,  strips  should  be  covered  in  situ  with  a  constant 
thickness  layer  of  high  thermal  conductivity  epoxy  that  will  also  water  proof  the  sensors, 
as  discussed  below. 

The  requirements  of  the  coating  and/or  coatings  upon  the  sensing  surfaces  are  severe. 
The  coating  must  have  a  high  dielectric  strength,  a  low  dielectric  constant  and  a  low  rate  of 
water  absorption  as  well  as  a  high  thermal  conductivity.  The  required  shear  strength  will 
be  relatively  low  for  an  adhesive,  thus  the  coating  will  only  need  a  nominal  value  for  shear 
strength.  The  coating  needs  to  be  applied  uniformily  over  the  sensor’s  surface  such  that 
directional  sensitivity  is  not  induced;  furthermore,  each  individual  sensor  needs  the  same 
coating  thickness  to  maintain  relatively  uniform  sensitivity  between  sensors.  Different 
coating  application  methods  were  investigated  with  cost,  uniformity,  repeatability  and 
relative  thermal  resistance  (tc/kc)  being  of  prime  concern. 

One  option  consists  of  using  a  drip  method  using  thinned  epoxies  with  high  thermal 
conductivities  (Master  Bond*  Polymer  Adhesive  EP  21  A,  25  BTU-in/(R-hr-ft^)).  The 
sensor  strip  is  mounted  vertically  with  a  layer  of  epoxy  applied  above  the  sensing  surface. 
Gravity  forces  the  epoxy  over  the  sensing  surface.  This  method  requires  the  strips  to  be 
mounted  to  the  model  before  the  coating  is  to  be  applied,  which  could  be  a  problem  due 
to  the  geometry  constraints.  Uniformity  could  be  easily  controlled  to  approximately  +/- 
20%  thickness  but  a  coating  of  no  less  than  3  mils  could  be  obtained.  The  thickness  of 
the  coating  could  be  reduced  if  extreme  care  was  taken  but  uniformity  between  sensors 
and  especially  between  sensor  strips  would  always  be  suspect.  This  method  should  be  only 
used  to  evaluate  different  coatings. 

In  the  second  method  investigated,  a  highly  thinned  layer  of  high  thermal  conductivity 
epoxy  was  sprayed  onto  the  sensing  surface  with  an  air  brush,  specially  made  to  shoot 
epoxy.  It  has  been  shown  that  the  thinnest  layer  which  may  be  applied  is  2  mils.  Better 
uniformity  is  achieved  with  this  method  over  the  drip  method.  However,  substantial  skill 
would  have  to  be  developed.  This  method  also  requires  the  bonding  of  the  strip  to  the 
model  before  the  coating  application,  but  does  not  require  the  strip  to  be  at  a  specific 
orientation.  Room  conditions  would  have  to  be  controlled  though. 

The  third  and  most  reliable  method  produces  a  constant  thickness  layer  of  high  k 
epoxy  without  special  application  skills  or  handling  of  sensors  before  their  mounting  on 
the  model.  A  strip  of  thin  (2  mils)  plastic  tape  is  applied  on  each  side  of  a  row  of  sensors 
after  they  have  been  mounted  on  the  model.  The  high  k  epoxy  is  applied  to  the  sensors. 
A  teflon  spreader  is  moved  across  the  top  of  both  plastic  tape  strips,  leaving  a  uniform 
layer  of  the  epoxy  in  between.  After  the  epoxy  partially  sets  up,  the  two  tape  strips  are 
removed,  along  with  any  epoxy  remaining  on  top  of  those  strips. 


*  Master  Bond,  P.O.  Box  522,  Teaneck,  N.J.  07666. 
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Another  possible  covering  is  a  conformal  polymer  coating.  The  best  conformal  coating 
found  for  this  application  has  been  Parylene,  applied  by  Nova  TVan  Corp.**  which  is  a 
subsidiary  of  Union  Carbide  Corp.  This  method  will  produce  uniform  coatings  which  will 
be  consistent  among  sensor  strips  since  the  group  of  strips  which  will  be  used  on  a  given 
model  will  be  coated  at  the  same  time  under  the  same  conditions.  A  coating  has  been 
investigated  which  when  applied  at  a  1  micron  thickness  will  give  a  breakdown  voltage  of 
greater  than  500  volts,  giving  an  order  of  magnitude  of  margin  upon  the  breakdown  voltage. 
The  coating  is  applied  such  that  no  pinholes  exist  at  a  0.1  micron  thickness.  Furthermore, 
the  water  absorption  rate  is  less  than  0.1%  by  mass  over  a  24-hour  period;  thus  substantial 
waterproofing  should  be  obtained.  The  thermal  conductivity  of  the  material  is  only  0.6 
BTU-in/(R-hr-ft^);  however,  due  to  the  thickness  of  application  a  reduction  of  an  order  of 
magnitude  of  effective  thermal  resistance  can  be  obtained  when  compared  with  a  sensor 
strip  manufactured  with  a  kapton  cover.  This  covering  must  be  applied  at  the  Nova  Trans 
plant  before  mounting  to  the  model.  The  conformal  coating  is  highly  flexible  and  will  not 
hinder  the  application  of  the  sensor  strips  to  the  model.  The  manufacturer  of  the  coating 
has  suggested  a  thickness  of  no  less  than  1  micron  such  that  the  coating  will  be  able  to 
sustain  mild  abuse.  The  constant  thickness  high  k  epoxy  method  was  still  judged  superior 
to  paralyne  since  in  situ  repairs  can  be  made  with  the  epoxy. 

The  epoxy  used  to  bond  the  sensor  strips  and  connectors  to  the  model  is  Hydrobond 
600,  which  is  manufactured  by  Superior  Polymer  Products***.  The  Hydrobond  has  a  water 
absorption  rate  of  0.1%  by  mass  over  a  24  hour  period  even  when  applied  underwater. 
This  epoxy  is  semi-flexible  when  cured.  The  user  is  able  to  control  the  pot  life  and  the 
viscosity  by  adding  and  subtracting  various  components.  The  hydrobond  is  also  used  to 
coat  the  connections  to  the  sensor  strip  since  lab  tests  have  shown  that  it  has  sufficient 
electrical  strength  for  applications  with  sensors  underwater.  The  electrical  properties  of 
the  connector  coating  are  less  stringent  than  the  sensor  coating  because  the  thickness  of 
the  coating  can  be  much  thicker  for  this  application.  Additional  electrical  insulation  can 
be  obtained  at  the  connections  by  using  strain  gage  technology,  which  incorporates  teflon 
sheets,  although  we  found  this  to  be  unnecessary. 

VI  Description  of  the  Constant  Current  Kilo  Anemometer  Electronics.  Data  Acquisition 
Hardware  fc  Program 

Our  purpose  is  to  apply  flush  mounted  hot  film  sensors  to  the  surface  of  a  model  to 
detect  the  location  of  the  three-dimensional  separation  lines.  Hot-film  sensors  are  a  good 
choice  because  they  are  directly  sensitive  to  the  local  shear  stress  and  relatively  easy  to 
manufacture  in  large  arrays.  For  the  application  to  the  16’  rotating  arm  submarine  model, 
on  the  order  of  1000  sensors  are  needed  to  define  the  position  of  a  separation  line  with 
sufficient  peripheral  spatial  resolution. 

The  most  practical  way  to  connect  1000+  constant  current  hot-film  sensors  is  in  series 
in  groups  of  N  sensors,  passing  a  common  constant  current  through  all  sensors  in  each 

**  Nova  Trans  Corp.,  100  Deposition  Drive,  Clear  Lake,  Wis.  54005. 

***  Superior  Polymer  Products,  1700  Duncan  Ave.,  Hubbell,  Mich.  49934. 


group.  This  reduces  the  number  of  leads  and  the  number  of  current  sources.  The  number 
of  leads  is  thus  reduced  to  N+2  instead  of  2N  for  each  group  (or  sensor  string).  Practical 
values  of  N  range  from  approximately  10  to  60,  depending  on  the  sensor  size  and  resistance, 
the  maximum  current  available,  the  maximum  voltage  available,  and  the  fluid  and  flow 
properties.  Such  systems  for  use  in  air  with  up  to  96  sensors  have  been  built  (Feyzi  et  al., 
1989). 

Here,  an  anemometer  system  for  up  to  1024  sensors  based  on  this  scheme  was  designed 
and  built  for  use  in  water.  The  electronics  of  the  system  consist  of  constant  current  sources, 
preamplifiers,  multiplexers,  a  computer  interface,  and  a  pair  of  transmission  line  interfaces. 
The  system  was  designed  to  fill  a  total  of  ten  circuit  boards,  eight  anemometer  boards  and 
two  interfacing  boards.  The  distribution  of  the  above  functional  sections  on  the  boards 
deserves  comment.  Since  multiple  printed  circuit  boards  were  required,  each  of  the  eight 
main  boards  was  designed  to  include  most  of  the  functions  for  128  of  the  sensors.  This  not 
only  allowed  a  variation  of  the  number  of  sensors  by  varying  the  number  of  boards,  but  it 
minimized  the  number  of  interboard  connections,  since  these  now  followed  the  first  stage 
of  16:1  multiplexers. 

The  etch  patterns  for  the  main  printed  circuit  boards  were  produced  at  Virginia  Tech. 
Prom  these  the  boards  were  manufactured  by  Appalachian  Electronics.  Only  one  board 
was  needed  for  the  140  foot  towing  tank  tests  and  the  VPI  and  SU  Water  Tunnel  tests 
described  below. 


A.  Constant  current  sources 

The  voltage-controlled  current  sources  consist  of  two  stages  of  op  amps  with  a  power 
PET  output  capable  of  producing  up  to  one  amp  of  regulated  current  at  up  to  40  v.  Each 
of  the  eight  boards  has  four  current  sources  (Walker  et  al.,  1992).  All  of  the  current  sources 
follow  the  same  calibration  curve  (within  one  per  cent).  Calibration  of  one  of  the  current 
sources  is  shown  in  Figure  7. 

The  current  is  software-controlled  via  the  voltage  from  a  D/A  output  of  the  computer. 
The  sensors  are  not  turned  on  until  the  model  is  in  motion.  The  current  can  then  be  fixed 
or  stepped  through  a  set  of  programmed  values  and  returned  to  zero  before  the  model 
is  brought  to  rest  and  the  test  terminated.  The  use  of  multiple  currents  for  each  sensor 
aids  the  correction  for  any  natural  convection  by  allowing  extrapolation  to  zero  current  and 
provides  a  least  squares  method  for  finding  the  forced  convection  heat  transfer  conductance. 

Each  current  source  operating  from  40  v  used  in  the  140  foot  towing  tank  provided 
power  for  up  to  14  sensors  at  450  ma  or  28  sensors  at  265  ma  with  (MINCO)  sensors.  The 
sensor  design  for  the  rotating  arm  test  will  allow  32  sensors  per  string  to  carry  375  ma. 


B. 


Anemometer  amplifiers 
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Each  anemometer  board  contains  128  differential  amplifiers,  consisting  of  32  PMI 
OP400  quad  op  amps  and  72  resistor  arrays  (Walker  et  al.,  1992).  Amplifiers,  such  as 
these,  are  required  to  eliminate  the  common-mode  voltage  resulting  from  the  other  sensors 
in  each  string.  Each  amplifier  observes  the  voltage  difference  across  one  of  the  sensors. 
The  amplifiers  also  provide  a  modest  gain  (3-10)  to  raise  the  signal  levels  away  from  noise 
and  interference. 

At  very  low  current  settings-  used  to  measure  the  sensor  cold  resistances,  typically  20 
ma-  the  current  load  due  to  the  amplifier  resistors  cannot  be  neglected.  To  correct  for  this 
an  algorithm  for  individual  sensor  currents  was  applied. 

C.  Analog  multiplexing/  digital  control 

Each  board  had  eight  16:1  multiplexers  (Harris  IH6116  for  the  towing  tank,  DG506A 
for  the  rotating  arm  system).  The  outputs  of  the  eight  multiplexers  from  each  of  the 
anemometer  boards  are  fed  through  four  more  16:1  multiplexers  on  a  ninth  board,  which 
reduces  the  1024  total  possible  sensor  signals  to  four.  The  four  remaining  signals  are 
transmitted  to  the  remote  computer  A/D  converter  board. 

To  control  which  sensor  is  read,  each  sensor  may  be  polled  independently  in  any 
order  and  at  any  rate  up  to  the  maximum  transfer  rate  by  the  host  computer.  Individual 
sensors  are  addressed  in  blocks  of  four  by  the  computer  via  an  eight  bit  address,  digitally 
transmitted  to  the  anemometer  boards.  Each  of  the  four  sensors  in  a  block  axe  distinguished 
by  its  A/D  input  channel.  For  the  present  application  all  of  the  sensors  are  polled  in  a 
fixed  sequence. 

The  maximum  rate  that  the  sensor  addresses  can  be  sequenced  is  much  higher  than 
the  rate  that  the  sensor  data  can  be  A/D  converted  and  stored  and  is  therefore  not  a  factor 
at  present. 

D.  Analog  transmission 

The  data  from  the  sensors  can  be  sampled  with  an  A/D  converter  inside  a  computer 
at  a  distance  of  250  feet  from  the  submarine  model,  although  in  the  140’  Basin  tests  the 
computer  was  mounted  on  the  carriage  with  no  more  than  20  feet  of  cable.  To  transfer 
data  at  the  required  rate  and  to  avoid  degradation  of  the  signal  due  to  interference  or  cable 
ringing,  an  analog  line  driver  was  designed.  The  transmission  lines  will  be  shielded  cables 
and  both  ends  of  the  transmission  lines  are  impedance  matched.  Tests  were  conducted  to 
verify  that  the  required  data  transfer  rate  of  10  samples  per  second  per  sensor  would  be 
reached.  The  tests,  conducted  over  both  a  coaxial  cable  and  shielded  twisted  pairs  nearly 
1000  feet  long,  demonstrated  that  data  rates  even  higher  than  the  required  rate  could  be 
achieved  over  the  required  transmission  distance. 

E.  Submersion  housing  for  electronics 

For  the  future  rotating  arm  tests  the  1000  anemometer  system  must  reside  inside  the 


# 
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flooded  maxine  model.  To  protect  the  electronic  boards  from  contact  with  water  or  the 
supporting  sting,  a  plexiglas  and  rubber  housing  has  been  designed,  machined  and  tested 
for  flt  in  the  model.  The  housing  consists  of  a  plexiglas  “box”  which  contains  a  neoprene 
liner  filled  with  an  insulating  gel  in  which  the  pc  boards  are  inserted.  The  gel  assures  that 
the  pressure  inside  and  outside  the  bag  is  the  same  at  any  submersion  depth  and  it  help® 
to  remove  heat  from  the  circuits.  The  housing  is  horse-shoe  shaped  so  it  can  wrap  around 
the  sting  and  to  allow  its  installation  from  the  side  of  the  model  while  the  sting  is  in  place. 

F.  Data  Acquisition  Hardware  and  Program 

Data  acquisition  was  performed  on  an  IBM  PS/2  Model  60  with  a  National  Instru¬ 
ments  MC-MlO-16  data  acquisition  board.  The  MC-MlO-16  contains  a  12  bit  analog-to- 
digital  (A/D)  converter  with  16  analog  inputs,  two  12  bit  digital-to-analog  (D/A)  convert¬ 
ers  with  voltage  outputs,  and  eight  digital  I/O  lines.  Multi-channel  data  acquisition  can  be 
performed  at  rates  up  to  100  KHz.  The  board  uses  computer  memory  as  a  storage  buffer. 

The  data  acquisition  program  is  written  in  C  language  to  maximize  data  acquisition 
rates  and  computer  storage.  Four  channels  on  the  MC-MlO-16  board  are  used  for  data 
acquisition,  the  eight  digital  lines  are  used  to  control  the  multiplexing  on  the  kilo-hotwire 
board,  and  the  voltage  output  is  used  to  control  the  current  through  the  sensors.  Currently 
the  data  acquisition  program  uses  buffered  I/O  with  a  maximum  data  size  of  128  Kbytes 
or  64  K  samples  of  data.  The  data  are  stored  in  binary  data  files  for  later  processing. 

The  data  reduction  program  is  also  written  in  the  C  language.  The  data  are  averaged 
for  each  sensor  over  the  total  number  of  data  points  taken;  Chauvenet’s  criterion  (Holman, 
1989)  can  be  used  to  remove  spurious  data  points.  The  data  are  then  corrected  to  take  into 
account  the  slightly  different  current  in  each  sensor  due  to  the  loading  of  the  resistor  arrays 
and  to  also  correct  for  the  different  gains  and  offsets  in  each  of  the  amplifiers  associated 
with  each  sensor.  The  resulting  voltages  are  converted  to  sensor  foil  temperature  difference 
by  (T/  —  Too)  equation  (36)  and  heat  transfer  rate  q  by  equation  (37).  A  least  squares  fit  of 
q  vs.  {Tw  —  Too)  by  a  straight  line  for  the  several  low  currents  where  forced  convection  is 
dominant  yields  a  slope  that  is  proportional  to  the  forced  convection  conductance  h.  (See 
equations  (23)  and  (38).) 

FVom  these  results  the  location  of  minimum  forced  convection  heat  transfer  conduc¬ 
tance  h  is  located  by  using  a  quadratic  least  squares  fit  to  the  h  data  from  5  adjacent 
sensors.  This  locates  the  minimum  skin  friction  magnitude  since  tw  =  n{dU/dt)w  ~ 
according  to  equation  (14). 

VII.  Experimental  Results  of  the  Performance  of  the  Constant  Current  Anemometer 
System  in  Locating  Three-Dimensional  Separations  on  a  Circular  Cylinder  at  Angle  of 
Attack  in  the  VPIfcSU  Water  Tunnel 

A  series  of  tests  was  conducted  in  the  VPI&SU  water  tunnel  (2’  x  2’  crossection)  to 
examine  the  performance  of  the  constant  current  anemometer  (CCA)  system  for  locating 
3-D  separations.  The  photographs  in  Figure  8  show  these  tests  in  progress.  A  6.625  inches 
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diameter  circular  plastic  pipe  was  used  as  the  test  body  at  an  angle  of  attack  of  22  degrees. 
Two  rows  of  CCA  sensors  were  mounted  around  the  cylinder,  as  shown  in  the  photographs. 
The  approach  water  flow  was  at  2  feet  per  second.  Prom  the  work  of  Poll  (1985),  it  is  clear 
that  we  should  expect  a  primary  separation  near  105  degrees  peripheral  angle  measured 
from  the  windward  side,  with  turbulent  reattachment  near  115  degrees  peripheral  angle. 
A  second  separation,  this  time  turbulent,  should  occur  at  some  higher  peripheral  angle 
somewhere  near  125  to  130  degrees. 

At  the  same  time,  we  used  “mini-tufts”  that  were  glued  directly  to  the  cylinder  up¬ 
stream  of  the  rows  of  sensors  to  observe  the  boundary  layer  behavior  in  that  vicinity. 
These  “mini-tufts”  were  fluorescent-dyed  0.001  inches  diameter  filaments  that  were  glued 
at  one  end  to  the  cylinder  surface.  While  not  readily  seen  under  visible  light,  the  unsteady 
motion  of  these  small  filaments  can  be  observed  with  UV  illumination.  Other  researchers 
(Crowder,  1980)  claim  that  these  tufts  can  follow  unsteady  motions  up  to  30  Hz.  However, 
in  the  tests  described  here  these  tufts  remained  too  far  away  from  the  wall  to  detect  the 
convergence  of  the  flow  from  each  side  of  a  3-D  crossflow  separation.  As  shown  in  these 
photos  for  the  0  degrees  angle-of-attack  case,  these  mini-tufts  move  erratically  away  from 
the  wall  of  the  cylinder  and  could  not  pick  up  the  separation. 

Figure  9  shows  data  from  8  adjacent  CCA  sensors  for  the  heat  transfer  rate,  q,  vs.  the 
sensor  cover  temperature  difference  (Tu,  —  Too)  for  8  different  sensor  currents.  A  straight 
line  closely  fits  the  data  for  a  given  sensor,  except  for  very  low  currents.  The  slope  of  a  given 
straight  line  is  proportional  to  the  heat  transfer  convective  conductance,  h,  according  to 
equation  (38).  The  cube  of  h  is  proportional  to  the  local  magnitude  of  the  surface  shearing 
stress  and  is  less  than  ±  4%  uncertain.  Estimates  of  {dU/dy)w  from  equation  (14)  indicate 
that  free  convection  effects  are  negligible  since  R  in  equation  (27)  is  much  less  than  0.1. 

Figure  10  shows  a  plot  of  the  variation  of  the  magnitude  of  the  shearing  stress  around 
the  circular  cylinder  at  22  degrees  angle  of  attack.  Since  the  local  minimum  in  the  magni¬ 
tude  of  the  surface  shearing  stress  is  the  location  of  a  three-dimensional  separation,  figure 
11  shows  that  one  3-D  separation  occurs  near  105  degrees  peripheral  angle,  reattachment 
occurs  near  115  degrees,  and  a  second  separation  occurs  near  125  degrees.  Considering  the 
resolution  limits  of  these  tests  and  the  flow  blockage  effects  of  the  floor  and  walls  of  the 
water  tunnel,  the  agreement  with  the  results  of  Poll  are  quite  good. 


Vin  Experience  and  Results  of  the  Performance  of  the  Constant  Current  Anemometer 
System  in  Locating  Three-Dimensional  Separations  on  a  6: 1  Prolate  Spheroid  at  Angle  of 
Attack  in  the  David  Tavlor  Research  Center  140’  Towing  Basin. 

A  series  of  tests  was  conducted  in  the  DTRC  140’  Towing  Basin  (10’  wide  x  5’  water 
depth  X  140’  long)  to  examine  the  performance  of  the  constant  current  anemometer  (CCA) 
system  for  locating  3-D  separations  on  a  6:1  prolate  spheroid  model.  This  54”  long  model, 
shown  in  Figure  11,  is  a  duplicate  model  of  the  one  used  by  Ahn  and  Simpson  (1992) 
in  studying  the  location  of  three-dimensional  separations  in  the  VPI&SU  Stability  Wind 
Tunnel.  A  2.75”  diameter  sting  was  attached  inside  the  model.  Figure  12  shows  the  1/16” 
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diameter  trip  wire  located  at  xjC  =  0.2. 

Three  rows  of  CCA  MINCO  sensors  were  mounted  near  the  stern  of  the  leeside  of  the 
prolate  spheroid  model  at  xjC  of  0.7,  0.8,  and  0.9,  as  shown  in  the  photograph  in  Figure 
13.  Because  the  spacing  between  sensors  was  1/2”,  the  angular  spacing  between  the  center 
of  each  sensor  was  7.1°  at  the  first  location  (x/C  =  0.7)  and  10.0°  at  the  third  location 
(x/C  =  0.9).  As  mentioned  above.  Hydrobond  epoxy  was  used  to  mount  the  sensor  strips 
in  place.  Conventional  plastic  electrical  tape  was  strapped  tightly  around  the  edges  of 
each  strip  and  around  the  model,  providing  enough  elasticity  to  hold  the  strips  in  place 
while  curing  the  epoxy  with  conventional  heat  lamps.  This  epoxy  or  an  additional  layer 
of  “5-Minute”  epoxy  was  also  used  to  form  a  fillet  fairing  between  the  model  surface  and 
the  sensor  strip  surface.  The  taut  electrical  tape  provided  a  smooth  external  mold  for  this 
fillet. 

Ribbon  cables  were  soldered  to  the  lead  pads  and  then  epoxied  to  the  model  surface  as 
shown  in  Figures  11  and  13,  extending  toward  the  leeside  in  order  to  minimize  interference 
with  the  separation  processes  nearer  the  windward  side.  The  ribbon  cables  were  then 
taped  to  the  sting.  The  sting  was  then  screwed  onto  the  140’  Towing  Basin  carriage  model 
mount,  as  shown  in  Figure  14.  All  measurements  were  made  with  the  model  at  15°  angle 
of  attack  to  the  direction  of  the  carriage  motion.  The  ribbon  cable  was  taped  to  the  sting 
mount  and  connected  to  the  constant  current  source  and  amplifier  board  (Figure  15)  and 
computer  (Figure  16)  on  topside  of  the  carriage. 

During  each  run,  the  carriage  was  brought  from  rest  to  a  nearly  constant  speed  before 
data  were  acquired.  For  the  5  fps  runs,  it  took  about  4  seconds  to  reach  a  constant  speed, 
leaving  about  12  seconds  for  obtaining  data  before  braking  of  the  carriage.  For  the  7 
fps  runs,  it  took  about  5  seconds  to  reach  a  nearly  constant  speed  with  only  3.6  seconds 
remaining  for  taking  data.  Data  were  obtained  over  enough  constant  speed  runs  for  a  total 
record  length  of  at  least  12  seconds  for  each  sensor  current.  No  runs  were  possible  at  lower 
speeds  because  the  carriage  speed  control  failed  before  these  data  could  be  obtained. 

As  in  the  earlier  tests  with  the  circular  cylinder,  g  vs.  the  sensor  cover  temperature 
difference  (Tw  —  Too)  was  obtained  for  7  different  sensor  currents:  5,  200,  300,  350,  400, 
450,  and  500  milliamperes  (ma).  A  straight  line  closely  fits  the  data  at  a  given  current, 
except  for  the  lowest  current.  The  slope  of  a  given  straight  line  is  proportional  to  the  heat 
transfer  convective  conductance,  h.  The  cube  of  h  is  proportional  to  the  local  magnitude 
of  the  surface  shearing  stress  and  is  less  than  4%  uncertain. 

Figures  17-20  show  plots  of  the  variation  of  the  magnitude  of  the  shearing  stress 
around  the  prolate  spheroid  at  two  axial  locations  (x/C  =  0.7  and  0.9)  for  two  different 
carriage  speeds.  No  data  were  obtained  from  the  middle  row  of  sensors  (x/C  =  0.8)  which 
showed  no  sensitivity,  even  though  much  work  was  done  trying  to  find  the  source  of  the 
problem.  A  “Thompson’s  Waterseal”  coating  had  been  applied  previously  to  the  surface 
of  only  this  row  of  sensors,  leading  us  to  believe  that  the  effective  waterproofing  of  the 
original  Kapton  was  greatly  reduced  by  this  coating. 


Since  the  local  minimum  in  the  magnitude  of  the  surface  sheauring  stress  is  near  the 
location  of  a  three-dimensional  separation,  Figures  17  and  18  show  that  two  separations 
were  detected  on  the  leeside  (<f>  >  90°)  at  x/C  =  0.7,  while  an  apparent  laminar-to- 
turbulent  transition  was  detected  on  the  windward  side.  Although  Ahn  and  Simpson 
(1992)  did  not  report  separation  locations  for  the  same  Reynolds  numbers  as  obtained  in 
the  towing  basin.  Figures  17  and  18  still  show  good  agreement  of  the  two  sets  of  data, 
considering  the  limited  spatial  resolution  of  the  sensors. 

Figures  19  and  20  show  the  results  at  x/C  =  0.9.  Since  we  anticipated  leeside  separa¬ 
tions  near  90°  from  the  earlier  work  of  Ahn  (Figure  21),  we  had  mounted  available  sensors 
(Figure  13)  so  as  to  detect  the  primary  separation.  However,  Figure  19  shows  a  separation 
near  110°  and  an  apparent  laminar-to-turbulent  transition  near  50°  for  the  slower  case. 
Figure  20  also  shows  an  apparent  laminar-to-turbulent  transition  near  50°  for  the  higher 
speed  case,  but  the  shear  stress  magnitude  is  continuing  to  decrease  at  the  leeward  end 
of  the  sensor  strip.  Figure  A.l  in  the  Appendix  (albeit  30°  angle  of  attack  data)  shows 
that  near  the  stern  of  the  prolate  spheroid,  the  shear  stress  magnitude  greatly  decreases 
over  a  relatively  small  peripheral  angle  change  near  the  primary  separation.  This  feature 
is  repr-^duccd  in  Figure  19.  Figure  A.l  also  shows  that  the  shear  stress  magnitude  remains 
at  a  relatively  low  level  further  around  the  leeward  side  of  the  body. 

Nevertheless,  at  x/C  =  0.9  there  is  not  good  agreement  with  Ahn’s  data.  The  Reynolds 
number  differences  could  be  part  of  the  reason,  although  Ahn  and  Simpson  do  not  show 
much  effect  of  Reynolds  number  at  the  location.  The  data  of  Kreplin  et  al.  (1985)  for  15° 
angle  of  attack  at  a  Reynolds  number  of  6.53  x  10®  using  natural  transition  indicate  that 
shear  stress  minima  occur  about  15°  further  around  the  body  at  all  x/C  locations  than 
the  results  of  Ahn  and  Simpson. 

However,  the  most  likely  reason  for  this  apparent  discrepancy  with  the  oil  flow  work  of 
Ahn  and  Simpson  is  that  the  leads  and  ribbon  cable  for  the  last  row  of  sensors  were  placed 
on  the  windward  side  of  the  model.  Even  though  the  flow  passing  over  these  protuberances 
would  not  pass  over  sensors,  it  could  have  been  affected  by  them.  Also,  Ahn  and  Simpson 
used  a  2”  diameter  sting  through  the  stern  end  of  their  identical  model,  whereas  we  used 
a  2.75”  diameter  sting  in  the  currently  reported  tests. 

Although  all  of  the  above  discussed  data  were  obtained  after  the  sensors  and  the  flow 
speed  had  reached  steady  state  conditions,  the  transient  responses  of  these  sensors  is  of 
much  interest  because  of  their  intended  future  use  under  transient  conditions.  Tests  were 
conducted  to  determine  the  worst  case  response  times  of  the  hot-film  sensors.  Since  the 
current  was  computer  controlled,  it  was  possible  to  turn  on  the  current  at  time  zero  and 
then  sample  the  history  of  the  response  voltage  as  the  sensors  heated  up.  The  results  for 
the  MINCO  sensors  on  the  prolate  spheroid  model  at  zero  flow  indicate  that  there  are 
three  time  constants,  0.01  sec.,  0.121  sec.,  and  4.36  sec.,  instead  of  one  as  modeled  in 
the  analysis  beginning  with  equation  (41).  In  other  words,  there  are  three  exponentially 
decaying  terms  instead  of  just  one  as  given  in  the  analysis.  Also  found  were  the  weighting 
functions  for  these  exponentials,  0.905,  0.035,  and  0.06,  which  indicate  the  fraction  of 
the  thermal  capacity  lag  effect  in  the  balco  foil,  the  kapton  cover,  and  the  kapton  and 


fiberglass  substrate,  respectively.  Examination  of  the  analysis  beginning  with  equation 
(41)  reveals  that  the  zero  flow  case  is  the  worst  case  in  regard  to  the  transient  response 
•  since  the  governing  equations  are  non-linear.  The  presence  of  some  mean  flow  convective 

heat  transfer  reduces  the  time  constant. 

A  more  detailed  analysis  of  the  transient  behavior  has  been  undertaken  and  will  be 
included  in  the  final  system  report  (Walker  et  al.,  1992).  Users  will  be  able  to  acquire 
0  transient  response  data  for  each  sensor.  This  will  allow  us  to  correct  transient  data  obtained 

in  the  future  rotating  arm  experiments  for  thermal  capacitance  lags. 


IX  Conclusions 

It  has  been  shown  in  this  report  that  the  local  minimum  in  skin  friction  magnitude 
is  approximately  the  location  of  a  three-dimensional  separation  on  the  leeward  side  of  a 
body  with  crossflow  separation.  If  the  relative  crossflow  momentum  of  the  two  boundary 
layers  that  converge  at  separation  is  much  different,  a  single  scheme  is  proposed  in  the 
Appendix  for  estimating  the  location  of  the  separation,  which  will  not  coincide  with  the 
local  minimum  of  skin-friction  magnitude. 

The  relationship  between  the  skin-friction  magnitude  and  the  heat  transfer  rate  from 
a  constant  current  anemometer  sensor  was  derived  and  shows  that  the  skin  friction  varies 
with  the  third  power  of  the  convective  conductance.  Thus,  the  convective  conductance  is  a 
sensitive  parameter  for  the  detection  of  local  minima  in  the  skin-friction  magnitude.  The 
transient  response  of  each  sensor  can  be  measured  and  used  to  correct  measured  data  for 
time  lags  during  transient  flow  conditions. 

These  sensors  can  be  waterproofed  by  several  means,  but  the  bvjst  appears  to  be  a 
uniform  coating  of  a  high  thermal  conductivity  epoxy,  which  will  also  protect  the  sensors 
from  abrasion. 

Two  series  of  measurements  were  performed  to  prove  this  concept.  One  set  was 
performed  in  the  VPI&SU  Aerospace  and  Ocean  Engineering  Water  Tunnel  with  a  circular 
cylinder  at  22°  angle  of  attack.  The  second  set  was  performed  using  a  6:1  prolate  spheroid 
of  15°  angle  of  attack  in  the  David  Taylor  Research  Center  140’  Towing  Basin.  In  both 
cases,  measured  separation  locations  agreed  with  those  obtained  in  other  experimental 
situations,  except  where  flow  conditions  and  interference  were  beUeved  to  be  different. 

Walker  et  al.  (1992)  will  report  complete  documentation  on  all  electronic  circuits, 
sensor  designs,  and  other  details  required  for  the  operation  of  the  kilo  constant  current 
anemometer  for  the  DTRC  rotating  arm  tests. 

The  success  of  this  type  of  sensor  for  determining  the  relative  convective  conductance 
in  water  leads  us  to  suggest  that  QUANTITATIVE  convective  conductance  and  skin  fric¬ 
tion  results  can  be  obtained  if  calibration  of  each  TYPE  of  sensor  is  done  with  respect  to 
known  values.  One  possible  calibration  setup  would  be  to  use  a  fully-developed  pipe  flow 
of  water. 
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APPENDDC  A 


THE  RELATIONSHIP  OF  SURFACE  PARAMETER  PARAMETERS 
TO  THE  LOCATION  OF  THREE-DIMENSIONAL  SEPARATION  LINES 

Along  a  three-dimensional  separation  line  on  the  surface  of  a  body  rises  a  stream 
surface,  across  which  no  flow  from  one  side  of  the  separation  line  can  pass  to  the  other  side. 
Fluid  moves  away  from  the  wall  with  a  relatively  large  normal-to-wall  velocity  component 
V  at  the  separation  stream  surface.  As  also  shown  in  Figure  1,  surface  skin  friction 
lines  converge  on  each  side  of  the  separation  line,  as  well  as  converge  toward  one  another 
(Tobak  and  Peake,  1982).  We  observe  that  low  velocity  near  wall  fluid  from  each  side  of 
the  separation  line  flows  toward  the  separation  Une,  so  a  local  minimum  in  skin  friction 
magnitude  must  occur  in  the  vicinity  of  the  separation  line.  Here  we  examine  the  governing 
equations  and  other  features  of  the  shear  stress  in  order  to  relate  the  behavior  of  surface 
parameters  to  the  location  of  separation. 

We  can  examine  the  continuity  and  momentum  equations  on  the  surface  for  an  im- 
compressible  flow.  From  continuity  we  obtain 
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At  the  surface  the  shear  stress  vector  is  given  by 
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where  the  continuity  equation  has  been  used  to  obtain  the  last  expression.  The  direction 
of  the  surface  shear  stress  vector,  and  therefore  the  direction  of  the  shear  stress  vector 
coincident  with  the  separation  line,  is  given  by  the  expression 


# 


^  =  tan 

Derivatives  of  equation  (A.  2)  and  (A.3)  produce 


(A.6) 


M 


d^v 


dy^ 


(A.7) 


w 


Thus,  the  near  wall  normal-to-wall  velocity  component  can  be  expressed  in  a  Taylor  series 
using  only  surface  parameters  for  the  first  few  terms. 


V  = 


df^v 


dy^ 


y^  d^V 

2  dy3 


(A.8) 


or 


i-r  " 

2ti 


V  =  -s^(V-f„)y=  +  + 


6/i' 


(A.9) 


Clearly,  we  see  that  the  surface  shearing  stress  contributes  lower  order  information  to  the 
location  of  a  separation  line.  Equation  (A.5)  relates  the  second  derivative  of  the  velocity 
component  nomoal  to  the  surface  to  the  divergence  of  the  wall  shear  stress  obtained  on 
the  wall.  Only  the  third  and  higher  derivatives  of  the  wall  pressure  field  are  related  to 
the  velocity  component  normal  to  the  wall.  Therefore,  wall  pressure  measurements  are 
relatively  insensitive  to  the  phenomena  of  separation.  Only  “massive”  separations  are 
detected  from  wall  pressure  measurements  (Vollmers  et  al.,  1982). 


Note  that  even  with  the  wall  shear  stress,  one  must  have  complete  data  over  the 
surface  in  order  to  perform  the  divergence  operation  and  locate  where  the  normal-to-the- 
wall  velocity  component  is  relatively  large  and  thus  the  location  of  the  separation  line.  The 
impracticality  of  an  extremely  large  number  of  shear  stress  sensors  leads  us  to  consider  an 
APPROXIMATE  method  for  locating  sepeiration  lines  with  fewer  shear  stress  sensors. 


To  a  first  approximation,  the  separation  line  will  be  in  the  vicinity  of  the  local  mini¬ 
mum  in  the  shear  stress  magnitude.  We  cannot  rigorously  say  how  close  the  local  minimum 
is  to  the  actual  separation  line,  except  for  the  special  case  where  the  converging  shear  flows 
and  their  surface  shear  stress  distributions  on  each  side  of  the  minimum  are  symmetric. 
In  that  case  the  symmetry  requires  that  the  local  minimum  in  shear  stress  magnitude  be 
coincident  with  the  separation  line. 


An  example  will  serve  to  show  how  well  this  approximate  method  works  for  crossflow 
separation.  Figure  A.l  shows  distributions  of  the  surface  shearing  stress  magnitude  and 
direction  at  various  axial  locations  on  a  6.T  prolate  spheroid  at  30°  angle  of  attack  (Kreplin 
et  al.,  1982).  (One  must  not  confuse  the  local  minimum  that  occurs  near  transition  from 
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laminar  flow  to  turbulent  flow  on  the  windward  side  (0  <  90°)  as  a  three-dimensional 
separation.) 

Figure  A. 2  shows  the  locations  of  the  first  local  minima  in  surface  shear  stress  mag¬ 
nitude  on  the  leeward  side  as  a  function  of  axial  length  along  the  body  as  determined 
from  these  data  (Meier  and  Kreplin,  1980).  As  pointed  out  by  Meier  and  Kreplin,  if  one 
takes  into  account  that  the  oil  flow  pattern  must  lead  to  a  primary  flow  separation  at 
lower  circumferential  angles  0,  the  agreement  with  hot -film  shear  stress  results  is  good. 
However,  Meier  and  Kreplin  did  not  point  out  that  the  agreement  with  the  oil  flow  results 
was  best  where  the  shear  stress  magnitude  distributions  were  more  symmetric  about  the 
minimum,  i.e.,  for  0.3  <  a:/2a  <  0.7.  For  x/2a  >  0.8,  Figure  A. 2  shows  a  larger  deviation  of 
the  shear  stress  minimum  location  from  the  oil  flow  results  while  Figure  A.l  shows  rather 
unsymmetric  shear  stress  distributions  about  the  minimum. 

This  leads  us  to  propose  a  refinement  in  locating  the  separation.  Figure  A.3  illustrates 
two  cases,  the  one  where  the  shear  stress  distribution  is  symmetric  about  the  minimum, 
and  a  second  where  it  is  not.  In  the  latter  case,  it  is  clear  that  the  near  wall  velocities  of  the 
flow  converging  from  the  left  are  greater,  having  greater  momentum  toward  the  separation 
line  and  pushing  the  minimum  shear  stress  magnitude  toward  the  lower  velocity  right  side 
of  the  separation.  This  means  that  the  minimum  shear  stress  location  is  closer  to  the  lower 
velocity  fluid  than  the  actual  separation  line. 

It  appears  that  the  maximum  (5ru,/30)max  and  minimum  {dr^/dt^)  slopes  of  the  shear 
stress  magnitudes  nearest  the  minimum  in  shear  stress  magnitude  can  be  used  as  measures 
of  the  relative  momentum  of  the  two  converging  streams.  If  so,  then  we  can  propose  to 
locate  the  separation  line  where  the  maximum  and  minimum  slope  lines  converge,  as  shown 
in  Figure  A.3  and  given  by 

,  {twB-TWa)  +  Sb4>B  -  Sa  <t>A  /  A  ■,r>\ 

= - (Sfl  -  Sa) - 

Of  course,  for  cases  where  the  shear  stress  distributions  are  symmetric  about  the  minimum, 
such  as  case  1  in  Figure  A.3,  the  results  from  equation  (A.  10)  are  coincident  with  the 
minimum  in  shear  stress  magnitude.  Figure  A.2  shows  the  results  obtained  using  equation 
(A.  10)  for  the  Kreplin,  et  al.  data.  Note  the  better  agreement  with  the  oil  flow  results  for 
x/2a  >  0.6. 

In  conclusion,  to  determine  the  approximate  location  of  crossflow  separation,  one 
must  have  at  least  enough  shear  stress  sensors  around  the  perphery  of  a  model  to  define 
the  local  minimum  and  define  the  nearby  maximum  and  minimum  slopes  of  the  shear  stress 
magnitude. 
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Figure  A.  2 


Figure  A.3 


Comparison  of  separation  regions  obtained  from  oil  flow  patterns  and  sur¬ 
face  hot  film  measurements,  a  =  30°  (Meier  and  Kreplin,  1980).  Circles 
denote  location  of  local  minimum  in  skin  friction  magnitude.  Squares  denote 
separation  locations  calculated  by  equation  (A.  10)  and  the  method  shown 
schematically  in  Figure  A.3. 

Examples  of  variations  of  the  wall  shear  stress  magnitude  in  the  vicinity  of 
a  local  minimum: 

(a)  example  of  locally  symmetric  distribution  of  |tw|  around  a  shear  stress 
minimum. 

(b)  example  of  locally  asymmetric  distribution  of  \tw  |  around  a  shear  stress 
minimum. 


Figure  1.  Schematic  diagram  of  a  three-dimensional  open  separation  (Maskell,  1955) 


Experimental  measurements  of  the  magnitude  and  direction  of  the  mean  wall 
shear  stress  at  one  axial  location  on  a  6:1  prolate  spheroid  (Kreplin  et  al., 


"igure  3. 


Shape  of  the  balco  foil  sensing  element  used  in  the  MINCO  sensor  (8  times 
actual  size). 


Figure  4. 


Shai)e  of  the  Advance  Circuit  Systems  halco  foil  sensing  clement  (8  time: 
actual  size). 


Water  flow  over  sensor 


fiberglas  model 


Figure  5.  Schematic  of  the  sensor  cross-section. 


Table  1.  Some  properties  of  materials  used  in  the 
MINCO  sensors  as  supplied  by  MINCO. 


material 

resistivity, 

^n-cm 

thermal 

conductivity, 

W/m-°C 

temperature 
coefficient  of 
resistivity ,  °C'  ’ 

Kapton 

polyamide 

0.98 

30  to  60  X  10  '^ 

Balco  foil 

70%  Ni,30%  Fe 

19.9 

200. 

0.0051 

Pyralux 

adhesive 

10^°  to 

0.22 

<  10'^ 

Figure  6.  (a)  Strip  of  MINGO  sensors  (actual  size);  (b)  Strip  of  Advance  Circuit  Sys¬ 

tems  sensors  (actual  size). 


Calibration  of  Constant  Current  Source 
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Figure  7.  Current  source  calibration  used  in  the  DTRC  140’  Towing  Basin  tests. 


Figure  8a.  Photograph  of  ultra-violet  illuminated  “mini-tufts”  mounted  on  a  circular 
cylinder  at  0  degrees  angle-of-attack  in  the  VPI&SU  Water  Tunnel. 


Figure  8b.  Photograph  of  the  apparatus  in  the  VPI&SU  Water  Tunnel  for  the  circular 
cylinder  at  22  degrees  angle-of-attack  tests  of  the  constant  current  anemome¬ 
ter  separation  location  system. 


Photograph  of  the  6.625”  diameter  circular  cylinder  in  the  VPI&SU  Water 
Tunnel  at  22  degrees  angle  of  attack. 
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Experimental  data  for  the  circular  cylinder  at  22  degrees  angle  of  attack. 
Measured  heat  transfer  rate,  <7,  for  each  of  8  constant  current  ancmoincter 
sensors  is  plotted  versus  the  sensor  cover  temperature  difference,  {T,„  —  T,^  ), 
for  each  of  8  different  current  settings. 


Figure  10.  Experimental  results  of  the  relative  magnitude  of  the  surface  shearing  sIk's^; 

versus  the  peripheral  angle.  Note  the  location  of  a  separation  near  105°,  a 
reattachmcnt  near  115°,  and  a  second  separation  near  125°. 


Figure  11. 


Photograph  of  the  54”  long  6:1  prolate  spheroid  with  the  sting  shown  at  the 
right. 


Figu.fe  12. 


Photograph  of  the  nose  of  the  6:1  prolate  showing  the  1/16”  diameter  trip 
wire  located  at  x/C=0.2. 


Figure  13.  Photograph  of  three  strips  of  sensors  located  at  x/C=0.7,  0.8,  and  0.9  on  the 
6:1  prolate  spheroid. 


Figure  14.  Photograph  of  the  prolate  spheroid  model  mounted  on  the  sting  of  the  car¬ 
riage  for  the  DTRC  140’  Towing  Basin. 


Figure  15.  Photograph  of  the  constant  current  source  and  amplifier  board  used  in  the 
140’  Towing  Basin  tests. 


Figure  16.  Photograph  of  the  constant  current  anemometer  board  and  data  acquisition 
computer  mounted  on  the  carriage  of  the  140’  Towing  Basin. 
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Figure  17. 


Relative  magnitude  of  the  surface  shear  stress  at  x/C=0.7  at  peripheral  an¬ 
gles  measured  from  the  windward  side  for  the  5  fps  flow  over  the  prolate 
spheroid  at  15  angle  of  attack.  T  denotes  location  of  transition  on  the  wind¬ 
ward  side;  R\  denotes  reattachment  location;  S\  denotes  location  of  primary 
separation;  Sz  denotes  location  of  tertiary  separation.  Primary  and  tertiary 
separation  from  wind  tunnel  data  at  Re  =  2.85  x  10®  are  indicated  by  arrows 
and  angular  values  (Ahn  and  Simpson,  1992). 
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Figure  18.  Relative  magnitude  of  the  surface  shear  stress  at  x/C=  0.7  at  peripheral 
angles  measured  from  the  windward  side  for  the  7  fps  flow  over  the  prolate 
spheroid  at  15  angle  of  attack.  T  denotes  location  of  transition  on  the  wind¬ 
ward  side;  denotes  reattachment  location;  S\  denotes  location  of  primary 
separation;  S3  denotes  location  of  tertiary  separation.  Primary  and  tertiary 
separation  from  wind  tunnel  data  at  Re  =  2.85  x  10®  are  indicated  by  arrows 
and  angular  values  (Ahn  and  Simpsor,  1992). 
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Figiire  19.  Relative  magnitude  of  the  surcface  stress  at  x/C  =  0.9  at  peripheral  angles 
measured  from  the  windward  side  for  the  5  fps  flow  over  the  prolate  spheroid 
at  15  angle  of  attack.  T  denotes  location  of  transition  on  the  windward  side; 
R\  denotes  reattachment  location;  Si  denotes  location  of  primary  separation; 
S3  denotes  location  of  tertiary  separation.  Primary  and  tertiary  separation 
from  wind  tunnel  data  at  Re  =  2.85  x  10®  arc  indicated  by  arrows  and  angular 
values  (Aim  and  Simpson,  1992). 
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Figure  20.  Relative  magnitude  of  the  surcface  stress  at  x/C  =  0.9  at  peripheral  angles 
measured  from  the  windward  side  for  the  7  fps  flow  over  the  prolate  spheroid 
at  15  angle  of  attack.  T  denotrs  location  of  transition  on  the  windward  side; 
Ri  denotes  reattachment  location;  Si  denotes  location  of  primary  separation; 
S3  denotes  location  of  tertiary  separation.  Primary  and  tertiary  separation 
from  wind  tunnel  data  at  Re  =  2.85  x  lO^*  are  indicated  by  arrows  and  anguKar 
values  (Aim  and  Simpson,  1992). 
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Figure  A. 2  Comparison  of  separation  regions  obtained  from  oil  flow  patterns  and  sur¬ 

face  hot  film  measurements,  a  =  30°  (Meier  and  Kreplin,  1980).  Circles 
denote  location  of  local  minimum  in  skin  friction  magnitude.  Squares  denote 
separation  locations  calculated  by  equation  (A.  10)  and  the  method  shown 
schematically  in  Figure  A. 3. 


Figure  A. 3  Examples  of  variations  of  the  wall  shear  stress  magnitude  in  the  vicinity  of 
a  local  minimum: 

(a)  example  of  locally  symmetric  distribution  of  (riv|  around  a  shear  stress 
minimum. 

(b)  example  of  locally  asymmetric  distribution  of  |rvy|  around  a  shear  strr.ss 
minimum. 


